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SEMI-ANNUAL REPORT 


This progress report for grant NAG 5-145/Supplement 4 covers the period 
January through August, 1986. The reporting period includes July and August 
1986 as the result of a two-month no-cost extension, requested February 19, 
1986 and subsequently approved. Most of the research discussed herein 
concerns the pole tide, the oceanic response to the Chandler wobble. The 
reader should consult earlier progress reports for relevant background 
material. 

In my Supplement 3 research, I developed a self-consistent dynamic theory 
of the pole tide in non-global oceans (oceans with continents present). The 
major goal of Supplement 4 was to extend that theory to the case of turbulent 
oceans. By the beginning of this reporting period, a method had been found 
for expressing the tide current velocities in terms of the tide height — with 
all variables expanded in spherical harmonics; see previous reports for 
details. I was then able to combine all tide equations into a single, 
non-differential matrix equation Involving only the unknown tide height. This 
equation could be written 

+M*£; 

where T ■ (t”} is the collection of spherical harmonic coefficients of the 
tide height, and M p and M R are the prograde and retrograde amplitude 
components of the wobble that creates the tide; the matrix B and vectors b r , 
b R as found in the case of turbulent oceans are given in the appendix. 

The pole tide must be constrained so that no tidewater flows across 
continental boundaries. During the reporting period this constraint, 
rephrased as in previous Supplements' research, was derived for the case of 
turbulent oceans; with the tide velocities expressed in terms of the tide 
height, the result is 
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with IB , b^ , b R also given in the appendix. 

These two matrix equations were combined in the manner devised for the 
Supplement 3 equations ("method (iii)" discussed below); simple matrix 
inversion then yielded the constrained solution. Such analytical equations 
for a tide in turbulent, non-global oceans have never before been derived, and 
represent an important development in the field. Programs to construct and 
invert the matrix equations were written (by this time I had become somewhat 
proficient at vectorized FORTRAN) ; preliminary results have been obtained and 
will be discussed near the end of this report. 


Prior to this reporting period, my Supplement 3 research appeared to have 
demonstrated that the pole tide in non-global oceans (without turbulence) 
modifies the Chandler wobble severely, causing appreciable damping and also 
lengthening the Chandler period by about 50% more than a static 


("equilibrium") pole tide would. One of the two major difficulties with these 
results, however, was the failure of the tide solutions to globally conserve 
mass (see previous reports); equivalently, the degree 0, order 0 tide height 
coefficient (To) failed to vanish. Although such a problem is commonly 
encountered in most tide theories, I felt it important that my tide solutions 
have a much smaller magnitude T c . During the reporting period, three methods 
were developed to achieve that goal: 


(i) Out of the nominal 36 equations (tide + boundary condition equations) 
whose solution would yield Tj through degree and order 5, one of the equations 
was replaced with an additional constraint that X Tds - 0 (tide mass is 
conserved in the oceans), or equivalently 






where are the harmonic coefficients of the ocean function; 


(ii) Same as (i) but with a second equation replaced by the constraint 
J Tds-0 (total tide mass on land is zero), or equivalently, 
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(iii) The number of boundary condition equations included in the matrix 
equation was increased. 


The result of method (i) was to significantly reduce the magnitude of T°, 
changing the other t” to a lesser degree. Since the two extra constraints 
combined are equivalent to requiring T° = 0, it is perhaps not surprising (but 
also not very satisfying) that method (ii) achieved a zero T 0 . 


The most reliable results were obtained with method (iii). 36 tide 
equations and ^ 3 6 boundary condition equations, each involving the 36 unknowns 
(T° through Tg), were combined to produce an overdetermined matrix equation of 
the form 
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with B a 72 x 36 matrix (72 equations in 36 unknowns). A least-squares 
approich led to the solution 
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In addition to this method being physically sound and mathematically 
consistent, a much reduced T 0 was obtained. 


Method (iii) was chosen for the Supplement 4 tide solutions as well. 


The other major difficulty with the Supplement 3 results was a variety of 
characteristics of the results — insensitivity to imposed wobble frequency. 



"zonality" of tide height coefficients, and so on — that suggested something 
was "wrong" with the solution. This possibility was reinforced when I began, 
during this reporting period, to compute net tide heights associated with the 
coefficients; significantly non-zero tide heights were found on land. I thus 
devoted a full summer month to re-deriving all of the Supplement 3 theory; 
another month was spent checking the entire FORTRAN code. No errors were 
found in either the theory or the programming. By the end of the reporting 
period, I had begun to realize that most of the characteristics I had judged 
"undesirable" were in fact consistent with the tide being close to static! 

The resolution of these difficulties, which turned out to require only a 
reinterpretation of existing results plus a small, simple additional 
computation, will be discussed in the next semi-annual status report. 

Preliminary results of the Supplement 4 research showed the effects of 
turbulence to be small. The fact that the solutions shared the same 
characteristics as the Supplement 3 results — despite the theoretical 
differences in approach (see discussion above or see previous report) — 
strengthened my conviction of the need to simply reinterpret the Supplement 3 
results. 


During this reporting period, my grant-related activities also included 
supervision of an undergraduate senior, S. Webb. Her project related to my 
original grant topic, the rotation of the coupled ocean — solid earth system. 
Her goal, determining the length of day changes associated with the coupling 
strengths implied by my original work, was not achieved during the one 
semester of her efforts (plans for her to continue during Summer 1986 were 
cancelled when she got the opportunity to work at the Lunar and Planetary 
Sciences Institute in Houston). 

In May 1986 I attended the Spring American Geophysical Union meeting and 
presented a talk entitled "The Self-Consistent Dynamic Pole Tide in Non-global 
Oceans"; a copy of the abstract is attached. During summer 1986, my article, 
"New Aspects of the Equilibrium Pole Tide," coauthored with D. J. Steinberg, 
was published in GJRAS ; copies have already been sent to NASA. Final 
revisions were made during spring 1986 to the manuscript "Another Look at 
North Sea Pole Tide Dynamics" (coauthored with J. R. Preisig); it was then 
accepted for publication. 



The Self-Consistent Dynamic Pole Tide 
in Non-global Oceans 


S.R. DICKMAN (Dept, of Geological Sciences, State 

University of New York, Binghamton, NY 13901) 

The pole tide is the oceans’ response to the Chandler 
wobble. The possibility of a dynamic response was 
previously Investigated for global oceans, using 
Laplace tidal equations augmented with linearized 
bottom friction [Dickman 1983]. A spherical harmonic 
approach allowed the equations to be expressed as a 
single non-differential matrix equation, with the coef- 
ficients of the tide height as the unknowns. In order 
to obtain a self-consistent determination of the tide 
and its effects on wobble, this equation was solved 
simultaneously with the Llouville equation, which 
governs the modifications of the Chandler period and 
damping by the tide. 

The previous theory has now been extended to the case 
of non-global oceans. The matrix equation is compli- 
cated significantly by the variable ocean depth, which 
was also expanded in spherical harmonics. The equation 
remains linear in the wobble amplitude forcing the 
tide, however, and this argues against the possibility 
[Carter, 1981] of frequency modulation of the Chandler 
period by a dynamic oceanic response. 

With continents present, the tide must be constrained 
so that no water flows across coastlines; such con- 
straints can also be expressed in matrix form. 
Preliminary self-consistent solutions have been 
obtained for oceans described by very limited harmonic 
expansions, for a variety of bottom friction strengths. 
1) The solutions reduce correctly to those of Dickman 
1985 when the oceans are treated as global. 2) For 
some combinations of ocean basin configuration and 
friction, a rotational resonance intensifying wobble 
appears. 3) There are problems in globally conserving 
mass, though application of the boundary constraints 
lessens the problem. 4) The oceans appear to dissipate 
significant wobble energy, and lengthen the Chandler 
period by ~50% more than a static response would. 
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